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Abstract. By means of a nonlinear separation of the variables in the governing full set of the magnetohydrodynamic 
(MHD) equations for axisymmetric plasmas we analyse an exact model for magnetized and rotating outflows 
which are hotter and overpressured at their axis. These outflows start subsonically and subAlfvenically from 
the central gravitating source and its surrounding accretion disk. Subsequently, they accelerate thermally and 
magnetocentrifugally and thus cross the appropriate MHD critical points, reaching high values of the Alfven Mach 
number. Three types of solutions are found : (a) coUimated jet-type outflows from efficient magnetic rotators with 
the flow confined by the magnetic hoop stress; (b) radially expanding wind-type outflows analogous to the solar 
wind, from inefficient magnetic rotators or strongly overpressured sources; (c) terminated solutions with increasing 
amplitude of oscillations in the width of the beam. In contrast to previously studied underpressured outflows, the 
transition from coUimated jets to uncoUimated winds is not continuous in the appropriate parametric space with 
a gap where no stationary solution is found. Superfast at infinity solutions are filtered by three critical surfaces 
corresponding to the three known limiting characteristics or separatrices of MHD wind theory. CoUimated and 
terminated solutions cross the slow, Alfven and fast magneto-acoustic critical points. Radially expanding solutions 
cross the slow and Alfven critical points while the last boundary condition is imposed by requiring that the pressure 
vanishes at infinity. 
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1. Introduction 

A well known example which demonstrates analytically 
that astrophysical jets can be accelerated and coUi- 
mated magnetically is the Blandford & Payne (1982) 



model. This model has been shown to be actually the 
prototype of the wide family of the so-called radially 
self-similar disk wind-type outflows which has been 
recently reexamined analytically and numerically (e.g 
Ouyed fc Pudritz, 1997[ |Vlahakis fc Tsinganos, 1998 



Krasnopolsky et al., 1999[ |Casse fc Ferreira, 2000 



Usty ugova ct al., 2000| I Kra snopols ky et al., 2003 

l^oh et al., 2002,,Casse fc Keppens, 2004f! 



A complementary wide class of MHD outflow solu- 
tions, which quantitatively demonstrated the transition 
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of coUimated outflows from efflcient magnetic rotators 
to uncoUimated outflows from less efficient magnetic 
rotators, is self-similar in the meridional direction 
(see ISauty et al., 2002a| henceforth STT02, and ref- 
erences therein). This class of models may describe 
ordinary stellar winds, or coUimated outflows composed 
of a central jet core surrounded by a disk wind ( 
Tsinganos & Bogovalov, 2002). Although this model 
is somewhat similar in geometry to an X-wind (e.g. 
|Shu et al., 1994| |Shang et al., 2002| ), it nevertheless has 
some differences, such as that it consistently solves 
the full set of the MHD equations from the source to 
the far region and also that the connection between 
the disk and the magnetosphere is an X point rather 
than a fan of concentrated magnetic flux. This class of 
analytical models may also be compared to the corre- 
sponding relaxation states of recent numerical simulations 
(e.g. [KoTde et al., 1998| [Bogovalov fc Tsinganos , 1999, 
■Keppens fc Goedbloed, 2000| [Matt et al., 2003^ 

Koidc ct al., 2000, Koide, 2003| ) as is discussed in 
Sec.l6.1JI 
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In such meridionally self-similar models, one may ei- 
ther prescribe the poloidal structure of the streamlines, 
or assume a relationship between the radial and longitu- 
dinal components of the gas pressure gradient. The main 
properties of the first class of solutions which are asymp- 
totically coUimated are outlined in Trussoni et al. (1997; 
henceforth TTS97) wherein the essential role of rotation in 
getting cylindrical coUimation has been demonstrated. On 
the other hand, if the two components of the pressure gra- 
dient are related, the meridional structure of the stream- 
lines is self-consistently deduced from the solution of the 
full set of the MHD equations. Such rotating and mag- 
netized outflows with a spherically symmetric structure 
of the gas pressure may be asymptotically superAlfvenic 
with radial or coUimated fieldlines, depending on the effi- 
ciency of the magnetic rotator ( |Sauty fc Tsinganos, 1994| 
henceforth ST94). 

In Sauty et al. (1999; henceforth STT99) we extended 
the results of ST94 by performing an asymptotic analysis 
of the meridionally self-similar solutions for a non spheri- 
cally symmetric structure of the pressure. It was pointed 
out there that a superAlfvenic outflow may encounter dif- 
ferent asymptotic conditions where it can be thermally 
or magnetically confined, and thermally or centrifugally 
supported. 

Current-carrying underpressured flows with a pressure 
increasing as we move away from the axis, were studied in 
STT02. They were found to be either thermally or mag- 
netically cylindrically coUimated around their axis, de- 
pending on whether the efficiency of the magnetic rotator 
prevails or not to the thermal confinement, respectively. 
They have been shown to be well suited to describe var- 
ious astrophysical winds and jets (see [Lima et al., 200l| 
IMeUani, 20"oTl [Sauty et al., 2003D . 

We complete here this work by studying overpressured 
outflows, i.e., with a pressure decreasing away from the 
system axis. Such outflows can only be coUimated via 
magnetic stresses if the magnetic rotator is sufficiently effi- 
cient. Otherwise the flow structure attains asymptotically 
a radial configuration. We present complete solutions that 
connect the base of the flow with its superAlfvenic regime. 
In particular we investigate if, and under which condi- 
tions, the basal region can be matched to the asymptotic 
solutions outlined in STT99. Conversely to the previous 
study reported in STT02, the present analysis requires 
a very careful topological study of the MHD self-similar 
equations because of the presence of a second X-type mag- 
netosonic critical point. 

In the following section 2 and in order to establish the 
used notation we briefly review the assumptions, parame- 
ters, variables and mathematical structure of the present 
model. In section 3 the asymptotic behaviour of the so- 
lutions presented in STT99 is also very briefly outlined. 
The results and parametric study are presented in Section 
4 while the main properties of the three classes of cylin- 
drical, radial and terminated solutions are summarized in 
section 5. Finally, in section 6 we discuss the astrophysi- 



cal relevance of our results, in particular in relation to jets 
associated with young stellar objects. 

2. Governing equations for meridional self-similar 
outflows 

We summarize here the main assumptions of our merid- 
ionally (6*—) self-similar treatment of the MHD equations. 
More details can be found in STT94, STT99 and STT02. 

2.1. Summary of the basic assumptions 

The basic equations governing plasma outflows in the 
framework of ideal MHD are the momentum, mass and 
magnetic flux conservation equations, together with the 
frozen-in law for infinite conductivity and the first law 
of thermodynamics. First, with axisymmetry the poloidal 
component of the magnetic field can be derived from a 
magnetic fiux function A{r, 9) in spherical coordinates 



B 



VA 



X if . 



(1) 



Next, for steady flows, we judiciously specify the 
meridional dependences of the velocity and magnetic 
fields, as well as of the density and pressure, V, B, P 
and p respectively, which may be written as follows (for 
details see ST94, STT99 and STT02): 
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P{R,a) = |p,K'n(i?)[i + H + Po . 



(2) 
(3) 
(4) 
(5) 
(6) 

(7) 
(8) 
(9) 



In the above definitions, k, A, d and Pq are model pa- 
rameters while we have also used the dimensionless mag- 
netic flux function a{R,0) = 2 A{r,e)/rlB^. Note that 
conversely to STT99 and papers before, we use a more 
flexible definition of the pressure function, as in STT02. 
The dimensionless pressure along the polar axis is defined 



C. Sauty et al.: Nonradial and nonpolytropic astrophysical outflows VI. 




b) 




Fig. 1. In a) is shown a plot of the dimcnsionlcss asymptotic radius of the jet Goo and in b) of the wavelength of the 
oscillations Aosc. in units of r*/A, vs. e/2A^ for a representative value of the terminal pressure, IIoo — O.OI. Each curve 
is drawn for a constant value of k/2A^ which ranges from —0.1 to — 10~^. To the right of the dashed line is the domain 
of pressure supported and magnetically confined jets (/b^ + fvp) while to the left of the dashed line is the domain of 
magnetocentrifugal jets (/s^ + fv^) (STT99). 



as n(i?) within some free additive constant Pq. It does not 
appear in the final dynamical equations as they depend 
only on the pressure gradient. However, this constant can 
be adapted to the boundary conditions. In radial solutions 
where 11 goes asymptotically to zero, Po should also van- 
ish to ensure that the temperature does not diverge as the 
mass density also goes asymptotically to zero. Conversely, 
in cylindrically coUimated flows where 11 can be a nega- 
tive function, Po should be adjusted such that the total 
pressure remains everywhere positive. 

The square of the poloidal Alfven number 



p 



(10) 



is assumed to be solely a function of the radial distance. 
For convenience we have normalized all quantities at the 
Alfven surface along the rotation axis, r ~ r^. The di- 
mensionless radial distance is denoted by _R = r/r,, , while 
-B*, T4 and p* are the poloidal magnetic field, velocity and 
density along the polar axis at the Alfven radius r,, with 
^2 = B^/Attp^. 

A second assumption is that the cylindrical distance 
vd[R, a) of a poloidal fieldline from the axis is separable 
in the variables R and a, as zu'^ = G^{R)a, where G^{R) is 
the cross sectional area of a flux tube perpendicular to the 
symmetry axis, in units of the corresponding area at the 
Alfven distance. Then, the dimcnsionlcss magnetic flux 
function a{R, 9) is related to G{R) through the following 
expression 



r:' 



G^{R) 



sm 



(11) 



Finally, for homogeneity with the notations in ST94, 
STT99 and STT02, we have also introduced the function 
F{R), which is the negative logarithmic derivative of the 



well known expansion factor used in solar wind theory 
| |Kopp fc Holzer, 1976J : 



F{R) = 2 



1 



d In G(i?) 



dlni? 



(12) 



We recall that the value of F dcflnes the shape of the 
poloidal streamlines. For F{R) — the streamlines are 
radial, for F{R) > they are deflected towards the polar 
axis (with F — 2 corresponding to cylindrical collimation) 
while for F{R) < they flare towards the equatorial plane. 

2.2. Parameters and variables 

The model is controlled by the following four parameters. 

— The parameter S which governs the non spherically 
symmetric distribution of the density with a linear in- 
crease (or decrease) of the density when receding from 
the rotational axis for 5 > (5 < 0). 

— The parameter A which is related to the rotation of the 
poloidal streamlines at the Alfven surface R ~ 1. 

— The parameter k which controls the non spherically 
symmetric distribution of the pressure. For k < (k > 
0) the gas pressure decreases (increases) by moving 
away from the polar axis. In this paper we confine our 
attention to overpressured jets, i.e., when k < 0. 

— The gravitational field is written as 



GM 



IK^ 



i?2 



(13) 



where Ai is the central gravitating mass. As a conse- 
quence there is an extra parameter v which is the ratio 
of the escape and flow speeds at the Alfven surface on 
the polar axis {R — 1), 



GM 



(14) 
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With the assumed axisymmetry, the original system 
of the MHD equations reduces to two coupled partial dif- 
ferential equations for the density and the magnetic flux. 
Furthermore, with the self-similarity assmnption, the com- 
ponents of V and B can be written as functions of 6 and 
three functions of R, namely G{R), F{R) and M{R). In 
this way, the momentum conservation law reduces to three 
ordinary differential equations which together with Eq. 
((T^ can be solved for the four variables A'P{R), F{R), 
n(i?) and G{R) (see Appendix A). 

2.3. Efficiency of the magnetic rotator 

By integrating the momentum equation along a fieldline 
we obtain the conserved total energy flux density per unit 
of mass flux density. This is equal to the sum of the kinetic 
and gravitational energies, together with the enthalpy and 
net heating along a specific streamline. In the framework 
of the present meridionally self-similar model, the varia- 
tion of the energy across poloidal fieldlines gives an im- 
portant extra parameter (STT99): 
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(15) 



which is a constant on all streamlines (ST94) . 

Physically, e is related to the variation across the field- 
lines of the specific energy which is left available to coUi- 
mate the outflow once the thermal content converted into 
kinetic energy and into balancing gravity has been sub- 
tracted (STT99). 

We can express e/2X^ in terms of the conditions at the 
source boundary Tq (see STT99 for details). 



e 
2A2 



£^Poynt.,o + -£^R,o + A_Eq 



E\ 



MR 



(16) 



where -Emr is the energy of the magnetic rotator (see Eq. 
2.5a in STT99), -Epoynt.^o is the Poynting energy, _Er^ is 
the rotational energy at the base and Ai?Q is the excess 
or deficit on a nonpolar streamline compared to the polar 
one of the gravitational energy (per unit mass) which is 
not compensated by the thermal driving. 



Ai^a 



GM 
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To (a) 



To (pole) 



GM {5 - K)a 
To 1 + 5a 



(17) 



For e > coUimation is mainly provided by magnetic 
means, while for e < the outflow can be confined only 
by the thermal pressure gradient, something which is not 
possible for overpressured flows. Accordingly, in STT99 
we defined flows with positive or negative e as Efficient 
or Inefficient Magnetic Rotators, respectively (EMR or 
IMR). 

A solution is determined by the four parameters v, 
e, K and A. The parameter 6 can be deduced from the 
constraint imposed by the integral e, Eq. H15|) . which has 



the following expression at the Alfvenic singular surface 

iR = iy. 

e = (k _ Sy + \^{t^ + 1) - (1 - k) + F2/4 , (18) 

where t[= (2 — F^)/p] is given by Eq. (|A.8|I . 

3. Asymptotic behaviour of the solutions 

For R ^ 1 the asymptotic parameters of collimated out- 
flows {Foo — 2, Goo and Moo bounded) depend on the 
value of e. Force balance across the poloidal streamlines, 
fvp + Ib^ + fv^ = 0, with /vp, Ib^ and fv^ the pres- 
sure gradient, magnetic stress and centrifugal volumetric 
force, respectively, calculates Moo and Goo as functions of 
the parameters e/{2}?), k/{2X'^) and Hoo- 

The asymptotic properties of these self-similar winds 
have been discussed in detail in STT99, and here we briefly 
summarize their main features for the case of overpres- 
sured outflows (k < 0), some of which are displayed in 

Fig.m 

— Two main asymptotic regimes exist. In one the outflow 
is collimated by the pinching of the toroidal magnetic 
field (e > eiim > 0) and in the other the outflow ex- 
pands radially (e < eiim)- 

— For K —f we have collimation for any value of e > 0: 
there is no pressure gradient across the streamlines and 
the flow can be supported only by the centrifugal force 

Ub, + fv, = 0) 

— Magnetically collimated flows are supported either by 
the centrifugal force or by the thermal pressure. For 
a given set of values of k/2X^ and Hoo, solutions with 
an increasing asymptotic radius Goo are found to pass 
from centrifugally supported to pressure supported. 
For each regime and a given value of e/2X'^ two so- 
lutions are found with a different asymptotic radius 

Goo- 

— Collimated streamlines always show oscillations. 
This behaviour is consistent with the results 
found in more general, non self-similar treatments 
( |Vlahakis &: Tsinganos, 1998| ) and numerical simula- 
tions. However as we move along a given curve with 
a fixed value of k/{2X^) and in the direction of in- 
creasing Goo (see Fig.^, in the region of centrifugally 
supported flows, e/{2X^) decreases and then increases 
again. In the region of minimum €/(2X'^) where the 
two regimes of centrifugally supported solutions merge 
we see that the wavelength becomes imaginary. This 
suggests that in this region cylindrical asymptotics is 
unstable. 

In conclusion and within the present model, from the 
asymptotic analysis it turns out that overpressured merid- 
ionally self-similar outflows from IMR should always ex- 
pand radially with an asymptotically vanishing pressure 
Hoo = 0. Conversely, outflows from EMR should undergo 
a transition from radially expanding to cylindrically colli- 
mating, as the efficiency of the magnetic rotator increases. 
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Fig. 2. Typical example of a radially expanding solution. In a) Vr/Vif along the polar axis is plotted vs. R the 
radial distance in units of the polar Alfven radius. The dotted line corresponds to the slow critical point and the 
dot-dashed one to the Alfven point. In b) the shape of the poloidal streamlines is plotted for e/(2A^) = —0.05 and 
k/(2A^) = —0.0005, with the dotted line indicating a streamline which is not connected to the central star but to the 
surrounding accretion disk. The various volumetric forces acting along and perpendicular to a given streamline are 
also plotted vs. R in the two lower panels. In c) positive forces (upper half) acting along the flow are accelerating 
forces while negative ones (lower half) are decelerating ones. The inertial force is negative as it is the opposite of 
acceleration. In d) positive forces (upper half) tend to coUimate while negative forces (lower half) are decoUimating. 
Symbols of the forces are defined in the text. 



4. Numerical results 

4.1. Numerical technique 

As in STT02, using routines of the NAG scientific pack- 
age suitable for the treatment of stiff systems and the 
Runge-Kutta algorithm, Eqs. I|A.1|I - l|A.7|l and I^SJ are 
integrated upstream and downstream of the vicinity of the 
Alfven transition (i?i„ — l±d_R) with Mi„ = l±pd_R and 
G,„ = 1 ± (2 - F^n)dR {F,n ~ F^). The slope p of M at 
i? = 1 is given in Eq. IJA.8(I . We first integrate upstream 
tuning the value of Fin until we select the critical solution 
that smoothly crosses the singularity corresponding to the 
slow magnetosonic point and reaches the base of the wind 
Ro with M ^ 0. With this value of Fin we then integrate 
downstream to the asymptotic region (with R^o usually 
between 10'' and 10^). 



Then, if the solution tends to become asymptotically 
radial or paraboloidal with a non zero pressure, we find 
that the transverse gradient of the pressure dominates and 
forces the flow streamlines to eventually flare towards the 
equator (or the pole if pressure is negative) at a finite 
distance. Such a solution is terminated. As is well known 
from the Parker wind theory, a physically acceptable so- 
lution which obtains radial asymptotics should satisfy the 
correct boundary condition at infinity, namely that the 
pressure should go to zero there. For this reason, we sim- 
ply tune the value of the pressure Hin (~ H,) such that 
n(_R — !■ cxd) vanishes. The various forces acting along and 
normal to a poloidal streamline are indicated in Figs. 12 
13 21 as follows (see STT02): fg, for the gravitational volu- 
metric force, fv for the inertial volumetric force, fs for 
the poloidal magnetic volumetric force, and /vp , fv^ , fs^ 
as defined in section 3. 
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Fig. 3. Typical example of a terminated critical solution. The same plots to those of Fig. [5| are shown for e/2A^ — 0.05 
and hi/2}? = —0.00015. In a) the dotted line corresponds to the slow critical point, the dot-dashed one to the Alfven 
point and the dashed line to the fast critical point. Positive forces are accelerating in the upper part of c) and 
collimating in the upper part of d), while negative forces decelerate in c) and decollimate in d). 



If an extra critical point appears downstream of the 
Alfven transition, we tune the value of the pressure Ilm 
(ss n*) such that the solution crosses this second X-type 
critical point. After the critical point the solution is ei- 
ther coUimated or terminated. We can always adjust the 
value Po such that the total pressure remains positive ev- 
erywhere (e.g. Figs. El and ^. 

Finally, if the solution naturally collimates and does 
not cross the second critical point, we obtain two alter- 
natives. Either we choose Po = as in STT99 and tune 
the value of the pressure Hm (~ H,) such that 11 (i?) is 
positive everywhere. Or, we let Po 7^ as in STT02 (e.g. 
Figs.0}, and choose a value of Po > such that the total 
pressure remains positive everywhere. 

It is worth to note that for k < 0, even though the 
pressure along the polar axis is always positive, it becomes 
negative for those nonpolar streamlines which correspond 
to a > — 1/k, as it happened in TTS97. Thus, conversely 
to the K > collimated solutions analysed in STT02, 
the present solutions cannot be extended to all stream- 
lines away from the flow axis. This limitation however is 
expected since is well known that the meridionally self- 



similar solutions are more adapted to describe the flow 
close to its axis (cf. ST94) than far from it. 

4.2. Behaviour of the solutions with k and e 

In this subsection we have fixed v — 1 and A = 1 (with the 
value of 5 deduced from Eq. IH^ and analysed the trends 
of the solutions for different values of e/2A^, k/2A^ and 
n, . Even though this is a rather restricted set which does 
not exhaust the whole space of the parameters, it may 
nevertheless illustrate the main characteristics exhibited 
by the solutions. 

For a constant negative value of k and for e increasing 
from negative to positive values three types of solutions 
are successively found. In the k/2A^ vs. e/2A^ plane of 
Fig. O are sketched the three distinct asymptotic regimes 
of the outflow: 

— Radial asymptotics solutions. An example of such a so- 
lution is shown in Figs. (21 where we have plotted the 
dimcnsionlcss radial speed along the polar axis vs. the 
radial distance R, the shape of the poloidal streamlines 
and the volumetric forces acting along and perpendic- 
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Fig. 4. Typical example of a coUimated critical solution. The same plots to those of Fig. 13 for e/2A^ = 0.05 and 
k/2A^ = —0.000005. In a) the dotted line corresponds to the slow critical point (SMSS), the dot-dashed to the 
Alfven point and the dashed to the fast critical point (FMSS). Positive forces are accelerating in c) (upper part) and 
collimating in d) (upper part), while negative forces decelerate in c) and decoUimate in d). 
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Fig. 5. Regions of the radially expanding (full squares), terminated (crosses) and collimated solutions (stars) in the 
plane of e/2A^ and k/2A^. 
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ular to the flow. Such solutions correspond to the filled 
squares in the [e/2A^, k/2A^] plane displayed in Fig.[5l 

— Terminated oscillating solutions, with a typical exam- 
ple shown in Figs.O Terminated solutions correspond 
to the crosses of the [e/2A^, k/2A^] plane (Fig. EJ. 

— Cylindrical asymptotics solutions, with a typical ex- 
ample shown in Figs. 01 Asymptotically cylindrical so- 
lutions correspond to the stars of the [e/2A^, k/2A^] 
plane (Fig.EJ. 

The same evolutionary trend is found for a given pos- 
itive value of e by increasing k from negative values (see 
Fig.EJ- The character of the solutions changes from radial 
asymptotics to cylindrical ones through a transition where 
the solutions are terminated. 

If e < 0, we are in the regime of IMR and the trans- 
verse pressure gradient does not let the flow to cylindri- 
cally collimate. In fact, solutions with a nonvanishing pres- 
sure are either terminated with negative pressure, or they 
have excessive flaring {F = —2) such that all stream- 
lines close at the equator, a rather unphysical situation 
as discussed in ST94. A third type of solutions is found 
when the pressure at the Alfven surface 11^, is tuned such 
that Hoc vanishes, as we mentioned above. Then, as pre- 
dicted by the asymptotic analysis, the flow streamlines 
asymptotically expand radially, with the Alfven number 
increasing unboundedly far from the base while the flow 
speed is bounded (V; — > Va,Br — > R^,M — > R). For 
a given e, by decreasing k to more negative values (Fig. 
E}, the solutions with radial asymptotics have a decreas- 
ing terminal velocity and initial pressure 11^,. This can 
be understood as follows. The transverse pressure gra- 
dient is proportional to kII, and lower values of 11^, are 
needed to open the lines radially and at the same time as 
the flaring is higher the velocity is lower as discussed in 



Tsinganos & Sauty (1992) This last result is unexpected 



from polytropic wind theory where a larger flaring leads to 
larger velocities |Kopp fc Holzer, 1976| but this is precisely 
what has been observed for the fast component of the so- 
lar wind during the minimum and the maximum of the 
last solar cycles (see |Wang, T995l|Wang &: Sheeley, "2003| ) . 

For e > and k lower than some threshold value ki, 
the same behaviour is observed. Solutions with an asymp- 
totically vanishing pressure are radial. Other solutions are 
either flaring with Foe = — 2 or refocalizing on the axis if 
the pressure becomes negative and then flaring again with 
F — —2 

When K reaches the value ki, we have the transition 
from squares to crosses in Fig.|Sl Now a second X-type crit- 
ical point emerges in the super Alfvenic regime for R > 1. 
Numerically, this second critical point appears at a flnite 
distance and is not coming from inflnity upstream, as one 
would expect. We notice further that, assuming |k| <C 1, 
for radially expanding solutions the quantity R^/G^ re- 
mains bounded, such that V does not reverse sign in this 
case, cf. Eq. (A. 4). Conversely if the flow tends to be col- 
limated or becomes radial very slowly the above quantity 
rapidly increases with R ior R :^ 1, leading to the ap- 



pearance of this new singularity [I? = in Eq. (A. 4)]. 
This explains the emergence of the second critical point 
in such cases. 

For a given range of k- values, [ki, K2], the two unphys- 
ical families which flare with F^o = — 2 separate at a finite 
distance. In principle the only possible solution would be 
the critical one that crosses this new X-type critical point. 
However, it turns out that this critical solution always 
shows downwind of the position of the critical point oscil- 
lations of increasing amplitude and eventually terminates 
in a loop at a finite distance. The termination position gets 
closer to the critical point by further increasing k to less 
negative values and/or increasing e. This type of solutions 
may be physically unacceptable since they do not extend 
up to infinity, unless they are terminated by a shock, with 
a positive pressure (we discuss this point in more detail in 
the next paragraph). In this case, this extra critical point 
seems to be the first of two or more fast critical transi- 



tions which seem to appear in the Weber & Davis (1967) 
1-D solution topologies and also in the 2-D analysis of 
Heyvaerts fc Norman (1989), Similarly to those examples, 
the solutions loop back rather sharply, returning upstream 
(Fig. 13. 

If the efficiency of the magnetic rotator increases fur- 
ther and/or k gets larger than K2, then a new fam- 
ily of cylindrically coUimated solutions enters the pic- 
ture. Topologically, they appear once the turning point 
of the terminated solutions reaches the X-type critical 
point. They have the typical properties of those from an 
EMR with K > (see STT02). The critical solution itself 
changes as it becomes also cylindrically coUimated. 

The various families of solutions which exist in that 
case are as follows. First, with rather high pressures we 
have solutions which still ffare with F^o = —2. Second, 
with rather low pressures we have the solutions which loop 
back. And finally, in between those two families of solu- 
tions, for a given intermediate range of 11*, we have the 
third family of cylindrically coUimated ones. 

One member of this third family of solutions crosses 
the X-type singularity which is still present at the bor- 
der between cylindrical and looping solutions. The other 
members of this third family of solutions are noncritical. 
The critical solution is analogous to the limiting solution 
of STT02 for k > 0. 

Furthermore we know from the asymptotic analysis 
that different branches of solutions are present, corre- 
sponding to centrifugally or pressure supported ffows. The 
present numerical results show that only the configuration 
with the smallest transversal radius can be attained by the 
jet, which is supported by the centrifugal force. Thus, as 
predicted by the asymptotic analysis performed in STT99, 
all cylindrically coUimated solutions including the criti- 
cal ones have almost the same asymptotic behaviour. In 
other words the pressure plays a minor role in achieving 
the asymptotic configuration of these solutions. 
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5. Properties of the critical solutions 

5.1. Cylindrically colli mated solutions 

An interesting novel feature in the present cylindrically 
coUimated solutions is the appearance of two X-type crit- 
ical points within the flow domain, in addition to the 
Alfven critical point. The only other known case where 
a unique steady MHD outflow solution is filtered by three 
critical points is the case of a radially self-similar solution 
(Vlahakis et al., 2000, |Ferreira fc Casse, 2004| ). In general, 
at such critical points the bulk flow speed equals to one 
of the characteristic speeds in the problem. Hence, it is 
of physical interest to associate the flow speeds at these 
critical X-type points to some characteristic MHD speeds. 
In that connection, we flrst note that the present solu- 
tions posses the symmetries of meridional self-similarity 
and axial symmetry. Thus, in spherical coordinates (r, 0, 
(f), the self-similarity direction is 9 and the axisymme- 
try direction is (p. Therefore, a wave that preserves those 
two symmetries should propagate along the f-direction 
in the meridional plane. First, the incompressible Alfven 
mode propagates along the magnetic field (B) with velo- 
city Va and in the direction f of the poloidal plane with 
a phase speed Va,r = B ■ f/y/Anp. Thus, at the Alfven 
point we should have M = 1. And second, the compress- 
ible slow/fast MHD modes propagate in the direction f 
with a phase speed Vx = Vsiow,r, or, Vx = Vfast.r- which 
satisfy the quartic 



Vx'-V^x(K' + Cf)+CfKf,. = 0. 



(19) 



Hence, when the above equation is satisfied the governing 
Eqs. IjA.ip . 1A.2|I and Eq. (|A.3|I have X-type singularities 
and Vx = Vp ■ f. 

On the other hand, it is well known that in the MHD 
flow system there exist two hyperbolic regimes wherein 
characteristics exist: the inner, which is bounded by the 
cusp and the slow magnetosonic surfaces and the outer 
extending downstream of the fast magnetosonic point. 
Within each of those two hyperbolic regimes, there exists 
one limiting characteristic or separatrix surface: the slow 
magneto-acoustic separatrix surface (SMSS) inside the in- 
ner hyperbolic regime and the fast magneto-acoustic sep- 
aratrix surface (FMSS) inside the outer hyperbolic regime 
(Bogovalov, 1994, Tsinganos et al., 1996). The true crit- 
ical points are precisely found on these two separatrices. 
For example, in the case presented in Figs. 4, the SMSS 
is at i? = 0.751 while the FMSS is located at i? = 4150. 

In the underpressured solutions studied in STT02 we 
have found only the X-type critical point inside the in- 
ner hyperbolic regime wherein the radial outflow speed 
is Vr — Vsiotu.r (ST94). Now, in the regime of cylindri- 
cally coUimated solutions (domain with stars in Fig. EJ , 
there exists a unique critical solution that also crosses the 
second critical point. This solution always has negative 
values of H(i?) asymptotically. Thus, in order to have pos- 
itive values of the total pressure everywhere in the flow, 
we have to adjust Pq to some positive value (Eq. EJ. The 



closer to zero is k the more negative becomes the function 
H and the larger is the minimum value needed for Pq. 

The fast magnetosonic nature of this second critical 
point can be analysed by drawing the characteristics in the 
vicinity of this separatrix critical surface, provided that 
we are able to define there the sound speed. Alhough the 
sound speed is ill-defined in our model, we can nevertheless 
deduce its value at the critical surface by following the 



steps presented in|Tsinganos et al. (1996) 


^'- dp 


(R,a) 


V^dIl{AP,R) 
2 dAP 


M^;+7+i/( 

{R,a) l + 6a 



(20) 



(21) 



with H{R, a) a function which becomes zero at the critical 
surface. The condition at the critical point is equivalent 
to setting equal to zero the denominator [Ea. lA.4| of Eqs. 
(TOll and l|X3|) . 



{M^-l)[l + ^§^ 



4 



{l-G 



2\2 



(1-M2)2 



= 0. (22) 



Together with the previous definition of the sound speed 
and Eq. (jA.l|l . this can then be put into the form of Eq. 
(O with Vx = Vr. 

Moreover assuming that H{R) = in all space we 
can calculate the slopes of the two characteritics in the 
regions where the equations are hyperbolic. The results 
are displayed on Figs. El for the cylindrically coUimated 
solution of Figs. 01 We show that there is an inner hy- 
perbolic domain bounded downstream by the slow mag- 
netosonic transition as well as an outer hyperbolic domain 
bounded upstream by the fast magnetosonic surface. As 
in Tsinganos et al. (1996) the transition from hyperbolic 



to elliptic in Fig. |B^is very close to the critical transition 
SMSS such that it is impossible to distinguish between 
the two without zooming closer than what is done in this 
figure. Instead in Fig.Oj, the FMSS is clearly distinguish- 
able from the FMS. Thus, in both cases the magneto- 
acoustic separatrices clearly differ from the corresponding 
magneto-acoustic transitions in the same way the ergo- 
sphere differs from the event horizon of a rotating black 
hole (see jSauty et al., 2002b| ). 

5.2. Terminated solutions 

We have seen that the main feature of the terminated so- 
lutions is the onset, downstream of the FMSS, of oscilla- 
tions with growing amplitude. This is likely to be related 
to some instability that does affect the configuration of 
the outflow which thus cannot attain any steady configu- 
ration. This is confirmed by the perturbative analysis in 
STT99: for this set of parameters the asymptotic solutions 
show a turning point (see Fig. Q] left panel, and Figs. 1 
in STT99) and in this region the typical oscillations of 
the asymptotic streamlines are not present, namely their 
wavelengths become imaginary (see Fig. ^ right panel, 
and Fig. 8, Eqs. 5.4 - 5.5 in STT99). 
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Fig. 6. Slopes of the two families of characteristics of the coUimated critical solution of Figs. 0]in each of the two 
hyperboHc regimes of the problem. In a) the slow magnetoacoustic separatrix surface (SMSS) is at i?=0.751 just before 
the slow magnetoacoustic surface (SMS). In b) the fast magnetoacoustic separatrix surface (FMSS) is at i?=4158 above 
the fast magnetoacoustic surface (FMS) at about R = 3000. Arrows indicate the direction of MHD signal propagation 
while two Mach cones above and below the FMSS are also shown. 



The cylindrical solutions with the second critical point 
are extended to infinity because there is some freedom on 
IIoo. For the terminated solutions, cylindrical asymptots 
are forbidden because the cylindrical regime is unstable 
(infinite wavelength of the oscillations). The only remain- 
ing solution is for the lines to become radial. However, by 
crossing the second critical point the value of 11 is fixed, so 
there is little chance that the solution becomes radial or 
paraboloidal because the pressure does not vanish asymp- 
totically. By plotting the forces across the lines, we see 
that in fact at the turning point the pressure gradient 
dominates. It is balanced only by the curvature force of 
the poloidal velocity which creates the turning point and 
the termination of the solution. 



regime. During that transition phase the outflow natu- 
rally passes from a stage of oscillations in its radius, Mach 
number and other physical parameters. These oscillations 
can be understood as the result of the interplay of the 
pinching magnetic tension force and the resulting reac- 
tion by the flow (conservation of angular momentum). In 
consistency with the previous explanation, a conspicuous 
feature appearing in the asymptotically radially solutions 
is the lack of any oscillation of the poloidal streamlines, a 
result confirming the physical origin of these oscillations. 



5.3. Radial solutions 

In outflows which are launched dominantly thermally (cf. 
Fig. 2c), the streamlines start with a basically radial shape 
close to their base. Further away, if the outflow is strongly 
overpressured (k < ki), the asymptotics remains radial. 
Hence, there is no drastic change in the geometry of the 
streamlines. On the other hand, in mildly overpressured 
onflows (k > K2), the asymptotics changes to cylindri- 
cal (in the intermediate case ki < k < K2 the oscilla- 
tions are so strong that the solution terminates). Thus, 
for K > K2, there exists a transition region between the 
base and the asymptotic regime wherein the cylindrical 
geometry is finally obtained after the basal roughly radial 



In that connection, theoretical arguments and various 
analytical self-similar solutions have shown that a notable 
common feature of all self-consistent, self-similar MHD 
solutions which become finally cylindrically collimated is 
that the outflow passes from a stage of oscillations in its 
physical parameters ( |Vlahakis fc Tsinganos, 1997| ). Such 
oscillatory behavior of collimated outflows is not restricted 
to the few specific models examined so far, but instead it 
seems to be a rather general physical property of an MHD 
outflow which starts noncylindrically before it reaches col- 
limation. Note that the same feature of oscillations has 
been also found in non-self-similar simulations of outflows 
which start radially before the magnetic tension converts 
them to a cylindrical shape ([Tsinganos et al., 2003|). 
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6. Conclusions 

6.1. Summary of results 

In this paper we continued the analysis of Meridional Self- 
Similar Models (MSSM, hereafter) by confining our at- 
tention to the study of outflows with a density increasing 
away from the axis and towards the surrounding stream- 
lines [cf. Eq. ((HJ with 5 > 0] and with a pressure decreasing 
from the axis [cf. Eq. with k < 0] . In such overpressured 
outflows with a central dip in the density distribution, the 
temperature is strongly peaked at the axis relatively to 
the surrounding regions, more than in the underpressured 
outflows studied in STT02. 

We have been able to construct solutions describing 
outflows starting subsonically and subAlfvenically from 
the central gravitating source and its surrounding accre- 
tion disk and, after crossing the MHD critical points, 
reaching high values of the Alfven Mach number. 

In terms of asymptotic profiles three broad types of 
solutions are found: 

— (a) CoUimated jet-type outflows from EMR where 
the outflow is confined by the magnetic hoop stress, 
provided that they are not too overpressured (k not 
too negative, i.e., k > K2). Among those solutions, a 
unique one crosses three critical points: the SMSS, the 
Alfven and the FMSS. We analysed those solutions and 
have shown that the separatrices indeed correspond to 
the three familiar separatrices in MHD wind theory. 
This class of critical solutions exhibits at large dis- 
tances oscillations between over- and under-pressured 
flow (Fig.Ej). 

— (b) Radially expanding wind-type outflows, analogous 
to the solar wind, for all IMR [cf. Eq. (fTK|l with 
e < 0] or strongly overpressured sources from EMR 
(e > 0, K < Ki). Those solutions cross the slow and 
the Alfven point and the initial pressure is fixed by 
the outer boundary condition that the terminal pres- 
sure should become zero. They do not show any inter- 
mediate topology of a third critical point. 

— (c) Terminated solutions. Such solutions cross again 
the three MHD separatrices. The onset of the in- 
creasing amplitude oscillations and the termination of 
the solutions can be understood because cylindrical 
asymptotics were shown to be unstable. As for disk 
wind solutions crossing all critical points, shown in 
Vlahakis et al. (2000)| such terminated solutions can 



support terminal shocks. This is the opposite of the 
solutions of refocalizing disk wind usually used in 
models (e.g. |Ferro-Fontan fc Gomez de Castro, 20031 ). 
Nevertheless, our terminated solutions are usually hav- 
ing very high temperature in the far regime and thus 
they are unphysical. 

Although from the asymptotic analysis presented in 
STT99 we could not exclude the existence of paraboloidal 
asymptotics, we could not find numerically any such solu- 
tion except in the limiting case e ^ k ~ (ST94). We con- 
jecture that these paraboloidal solutions could be found 



in principle in the transition region between radial and 
cylindrical asymptotics. However steady equilibrium is im- 
possible in this region and the solutions terminate before 
reaching their asymptotic regime. 

6.2. Astrophysical implications 

Possibly related to the previous discussion on the lack of 
paraboloidal solutions, we note that, conversely to under- 
pressured jets studied in STT02, the transition of coUi- 
mated jets to uncoUimated winds is not continuous in the 
parametric space showing a gap where stationnary solu- 
tions do not exist. We are tempted to conjecture, then, 
for some jets with strong transient events, such as violent 
outbursts, the following scenario: if the outflow configu- 
ration is at the interface of the regime with coUimated 
and non coUimated solutions, an outburst could be asso- 
ciated with the flip over between the two different classes 
of asymptotically coUimated and radial solutions. 

One of the major outstanding questions in astrophysi- 
cal jets research is how they are generated. Take for exam- 
ple the case of the closer and thus better resolved jets as- 
sociated with young stellar objects. In that connection one 
may say that through a combination of observations and 
numerical simulations we do know several details about 
the propagation of these jets in the parent cloud and their 
interaction with their environment but we know relatively 
fewer details about their generation at the "central en- 
gine" (cf. IHartigan, 2003||Ray fc Bacciot ti, 2003). At the 
same time various studies have shown that a high degree 
of coUimation is already achieved very close to the source, 
namely at 10 to 20 AU ( |Woitas et al., 2002| ). 

In that content, two wide classes of models are 
available today to study analytically the launch- 
ing and eventual coUimation of MHD outflows 
( [Vlahakis fc Tsinganos , 1998); first, the family of 
the so-called radially self-similar models (RSSM) which 
have as their prototype the Blandford fc Payne (1982) 
model and second the family of the so-called meridionally 
self-similar models (MSSM) which have as their proto- 
type the Sauty & Tsinganos (1994) model and which we 
have explored in this series of papers. A third class of 
models are the so-called X-wind models where mass loss 
originates at a fan of concentrated magnetic flux in the 
inner disk radius ( |Shu et al., 1994||Shang et al., 2002D . 

In an analytical MHD treatment of the problem of out- 
flow launching and subsequent coUimation, in the RSSM 
the driving force and coUimation mechanism are basi- 
cally magnetic. Among the limitations of the RSSM is 
that they are invalid close to the jet axis where they 
have singularities and are thus more appropriate to de- 
scribe disk-winds which coUimate within several AU from 
the star ( [Ferreira, 1997| henceforth F97). Also, after sev- 
eral Alfven radii when the streamlines reach their maxi- 
mum cylindrical radius, they slowly refocus towards the 
axis and the solutions terminate. Another difficulty of the 
cold plasma RSSM is that they predict rather large termi- 
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nal speeds and too low densities and ionisation fractions 
and do not accomodate some efficient heating mechanism 
which is needed in order to explain the observed emission 
| |Dougados, 2003| . However, Casse & Ferreira (2000) have 
shown that by including a hot corona above the disk it 
would help to increase the mass loss rate and thus the 
terminal densities. In any case, the RSSM deal consis- 
tently with the accretion ejection problem and have been 
recently used with some success to compare observed jet 
widths and coUimation scales in several T Tauri micro jets 
for which we currently have the corresponding observa- 
tions (Dougados et al., 2000, Pesenti et al., 2003). 

On the other hand, in the MSSM for jet acceleration 
and coUimation the driving force is a combination of ther- 
mal and magnetocentrifugal terms while coUimation can 
be also achieved by a combination of pressure gradients 
and magnetic tension forces. In the MSSM however, if the 
source region of the outUow is restricted to be only the 
stellar base, the resulting mass loss is unrealistically low, 
unless it includes the inner part of the disk. Nevertheless, 
an interesting fact is that observations clearly show that 
jets extend to relatively great distances (100 to 1000 AU) 
from the protostar where the observation of forbidden line 
emission means that the jets are still warm/hot at such 
large distances. But then if jets are launched from small 
regions and also expand, they should cool adiabatically. 
The question arises then on how do these jets remain hot 
at such distances from the protostar. Clearly a heating 
mechanism is needed. Hence, observations seem to suggest 
that thermal gradients, which may originate in a stellar or 
an accretion heated disk-corona, play an important role in 
accelerating the Uow (Dougados, 2003). Such a heating is 
a basic ingredient of the MSSM and may thus explain the 
puzzle. In addition steady MSSM stay tightly colUmated 
to unlimited distances from the source without a need to 
refocus towards their axis along which they are valid with- 
out any singularity. Furthermore, it has been suggested 
that HH jets may be the progenitors of the (uncoUimated) 
solar wind outflow, a form in which jets eventually evolve 
after the star looses angular momentum becoming an in- 
efficient magnetic rotator (ST94, STT99, STT02 and this 
paper). 

However, recent numerical simulations of magnetocen- 
trifugally coUimated outflows from a rotating central ob- 
ject and/or a Keplerian accretion disk have shown that 
relatively low mass and magnetic fluxes reside in the pro- 
duced jet as compared to the surrounding wind (Tsinganos 
& Bogovalov, 2002, IMatt et al., 20 03). This is also the 
case in the solutions presented in ST94 and STT02. In 
ST94 it is pointed that a signiflcant fraction of the to- 
tal mass loss rate of the jet is originating in the disk. 
Observations however indicate that in jets from young 
stellar objects, the colUmated outflow carries higher fluxes 
than these studies predict. As a solution to this problem 
it has been proposed that jets may be described as a two- 
component system composed of an outflow originating at a 
central object which is surrounded by a disk-wind (ST94, 
IKoide et al., 1998[ STT02, Tsinganos & Bogovalov, 2002). 



In that respect, Hartigan et al. (1995) have identifled an 
outer low velocity component (LVC) with velocities in the 
range of 10 to 50 km s^^ along with a high velocity compo- 
nent (HVC) with radial velocities of a few hundred km s^ ^ . 
According to Kwan & Tademaru (1995) the LVC is prob- 
ably a low-velocity disk-wind that encompasses the jet. 
This view was confirmed later by Bacciotti et al. (2002) 
using HST/STI data to investigate the velocity structure 
of the DC Tau jet. They found that the kinematics follows 
an "onion- like" structure with HVC closer to the jet axis 
and a LVC spread out wider. 

Double component jets are also clearly appearing in 
time dependent simulations of jets around black holes and 
stars (e.g. |Kudoh et al., 1998| |Koide et al., 1998| |. The 
magnetic field does not penetrate the black hole mag- 
netosphere, thus the inner plasma is compressed and a 
pressure driven outflow develops. The surrounding wind 
is centrifugally driven from the disk and magnetically col- 
Umated. However the main difference with double jet com- 
ponent simulations and the present analytical solutions is 
that the first describe dense core jets while the second 
model hollow jets. At this point note that both analytical 
(Hanasz et al., 2000| ) and numerical ( |Kudoh et al., 2002) 1 
studies of the flow stability tend to show that the inner 
flow is probably more stable than the Keplerian outer 
part and this is particularly true in the case of hollow 
jets. Further comparison of simulations and our analyt- 
ical modelling is however difficult because the boundary 
conditions generally used in disk winds are rather different 
from the self-similar assumptions used here. 

In accordance with the above theoretical and obser- 
vational difficulties encountered by the single-component 
models, we propose that jets may indeed be described as 
a two-component outflow system. The model presented in 
this paper and in STT02 is by itself a double component 
jet structure for cylindrically colUmated solutions. One 
part of the jet comes from the star itself. The other comes 
from the inner boundary of the disk which is connected 
with the stellar magnetosphere. Such double structure can 
be applied directly to model jets from T Tauri stars with 
low mass accretion rate, like RY Tau for instance. However 
for T Tauri stars having higher mass accretion rate, like 
DC Tau, a consistent model would be one with an inner 
outflow described by a ST94-type MSSM, surrounded by a 
wider F97-type RSSM disk-wind part. For example, in the 
present MSSM the velocity is peaked at the axis and the 
degree of coUimation increases with velocity, as observed. 

A supporting evidence for the previous scenarios comes 
from recent findings on the rotation of jets from T Tauri 
stars with high accretion rate, by using either near- 
infrared long-slit spectroscopy in a series of distant knots 
(several 1000 AU away), or, by using HST/STIS observa- 
tions mucU closer (within 100 AU) (Bacciotti et al., 2002, 
[Coffey et al., 20 03). The remarkable fact is that the mag- 
nitude of the only recently inferred toroidal velocity in 
the jets (5 - 15 km s~^ at distances of 20 to 30 AU from 
the flow axis and at around 100 AU from the plane of 
the disk) is precisely what some time ago was already in- 
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ferred from MHD jet launching models (see for example, 
|Tsinganos fc: Trussoni, 1991, Fig. |H| and Tsinganos et al. 
1992, Figs. 2,3 for a MSSM; or, Vlahakis et al., 2000, Figs. 
5,7 for a RSSM). It should be interesting then to further 
discuss new observations in the context of the two analyt- 
ical available models (MSSM and RSSM). 
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Appendix A: MHD Equations for meridionally 
self-similar flows 

Under the assumptions of axisymmetry and meridional 
self-similarity, the MHD equations reduce to the following 
three ordinary differential equations for n(i?), lvP{R) and 
F{R): 



dn 

dR 



2 
G4 



dR 






NPR^ 



dF{R) _ 7V>(i?, G, F, M2, U; k, S, v, A) 
dR " RV{R,G,F,M^;k,\) ' 

dA/2(i?) 7VA/(i?,G,i^,M2,n;K,(5,i^,A) 



dR 



RV{R,G,F,M^;k,\) 



(A.l) 
(A.2) 
(A.3) 
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where we have defined: 



X 



+ [2kUG^R^ + (F + 1){F - 2)] 



4^2 M^D^I' ^^-^^ 



+KllR^G^M^^ -^{F- 2){Ak% +F + A) 
2 4 G^ 

+ ^(F-2)(8«|J + F2 + 4^ + 8) 

-\^R\F-2)^ 

+A^i?^(2M^ + f-2)(§-^§). (A.6) 

with 

iVB = l-G^ A^y^M^-G^ i:> = l-M2. (A.7) 

The definitions of the various parameters is discussed in 
Sec. 2. 

At the Alfven radius, the slope of M^{R = 1) is 
p — (2 — F^,)/t, where r is a solution of the third degree 
polynomial: 

,3,o_2, r«n* , F^-A ^\_, {F,-2)F, 



and the star indicates values at i? = 1 (for details see 
ST94). 



